Thermal conductivity of supercooled water 
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The heat capacity of supercooled water, measured down to —37 °C, shows an anomalous increase 
as temperature decreases. The thermal diffusivity, i. e., the ratio of the thermal conductivity and the 
heat capacity per unit volume, shows a decrease. These anomalies may be associated with a hypoth- 
esized liquid-liquid critical point in supercooled water below the line of homogeneous nucleation. 
However, while the thermal conductivity is known to diverge at the vapor-liquid critical point due to 
critical density fluctuations, the thermal conductivity of supercooled water, calculated as the prod- 
uct of thermal diffusivity and heat capacity, does not show any sign of such an anomaly. We have 
used mode-coupling theory to investigate the possible effect of critical fluctuations on the thermal 
conductivity of supercooled water, and found that indeed any critical thermal-conductivity enhance- 
ment would be too small to be measurable at experimentally accessible temperatures. Moreover, the 
behavior of thermal conductivity can be explained by the observed anomalies of the thermodynamic 
properties. In particular, we show that thermal conductivity should go through a minimum when 
temperature is decreased, as Kumar and Stanley observed in the TIP5P model of water. We discuss 
physical reasons for the striking difference between the behavior of thermal conductivity in water 
near the vapor-liquid and liquid-liquid critical points. 



I. INTRODUCTION 

Supercooled water exhibits several thermodynamic 
anomalies, perhaps the best known of which is its den- 
sity maximum just above the freezing point, at 4 °C 
The isobaric heat capacity @, H|, isothermal compress- 
ibility Q, and thermal expansivity [l[ show anomalies 
that resemble critical-point power laws Moreover, 
the correlation length characterizing density fluctuations 
increases markedly upon supercooling 

In 1992, Poole et al. proposed a coherent and particu- 
larly fruitful account of the thermodynamic anomalies of 
water: that deep in the supercooled region is a first-order 
phase transition between two liquid states distinguished 
by their different densities. This transition line would 
terminate at a critical point, and proximity to this crit- 
ical point could explain the anomalous behavior of the 
response functions [6[ . The liquid-liquid phase transition 
and critical point are hypothesized to occur at temper- 
atures and pressures that are inaccessible to experiment 
due to unavoidable homogeneous nucleation of ice I h [l[ . 

This conjecture has given rise to several models, most 
of which propose two different ways in which water 
molecules might form hydrogen bonds. One of these 
super-molecular arrangements (the high-density liquid or 
HDL) is associated with higher density and is favored 
at higher temperatures and higher pressures; the other 
(the low-density liquid or LDL) is associated with lower 
density and is favored at lower temperatures and lower 
pressures. 

Intriguing experiments on the melting lines of 
metastable phases of D 2 and H 2 ice @, H[ have lent 
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additional support to this hypotheses, as has the find- 
ing of a phase transition in glassy water between two 
amorphous states of different density [|[ . A recent study 
comparing the locus of glass-transition temperatures as 
a function of pressure in simulations of water molecules 
and in real water has found the behavior of real water to 
bear a stronger resemblance to the model that exhibits a 
liauid-liquid phase transition [Io| ■ Other simulations [TiT - 
Il3| ] note an increase in tetrahedrality in water upon su- 
percooling, giving some idea as to the difference between 
the HDL and LDL states. Equations of state based on 
the hypothesized existence of a liquid-liquid critical point 
have provided increasingly accurate accounts of the ther- 
modynamic properties of water over the last several years 
[lj-[l6| . The most recent two-state model of Holten and 
Anisimov (with entropy-driven liquid-liquid separation) 
[l7| shows excellent agreement with the thermodynamic 
data with fewer adjustable parameters than any model so 
far. We have used the predictions of this two-state equa- 
tion of state (which we abbreviate "TSEOS" ) for ther- 
modynamic properties in our calculations of the thermal 
conductivity. 

While most of the phenomenology surrounding two- 
state models generally and the liquid-liquid critical point 
in particular has focused on thermodynamics, there are 
important implications for dynamics as well. To take one 
example, the viscosity of water decreases upon compres- 
sion, which is anomalous. A suggested explanation for 
this anomaly is that compression forces a greater fraction 
of water into the HDL state, which has greater mobility 
than the tetrahedrally ordered LDL state [l8|. Further- 
more, the dispersion of sound at high frequencies seems 
likely to reflect viscoelastic behavior associated with a 
structural relaxation in supercooled water [19(. 

In this paper we examine the thermal conductivity of 
supercooled water in light of the two-state conjecture, 
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and we examine the possible effects that critical fluctua- 
tions might have on thermal conductivity. According to 
our calculations, any effect of critical fluctuations asso- 
ciated with the virtual liquid-liquid critical point on the 
thermal conductivity would be too small to be measured 
at experimentally accessible temperatures. Remarkably, 
the behavior of thermal conductivity can be fully ex- 
plained by the anomalies of the thermodynamic proper- 
ties. The difference between the behavior of the thermal 
conductivity in the vapor-liquid and in the hypothesized 
liquid-liquid critical regions is the result of differences in 
both the dynamic and thermodynamic environments in 
the respective regions. 



II. THERMAL DIFFUSIVITY AND THERMAL 
CONDUCTIVITY 

The thermal diffusivity a of water has been measured 
by Taschin et al. down to 256 K [20j and by Benchikh et 
al. down to 250 K [2l[ , and it decreases steadily with de- 
creasing temperature, as shown in Fig. 1. Thermal con- 
ductivity A can be calculated from these data by means 
of the formula A = pc p a, since isobaric specific heat ca- 
pacity c p and mass density p are both known experi- 
mentally in the relevant temperature range. The heat 
capacity changes little in the range for which thermal- 
diffusivity data are available (Fig. 2), so in that temper- 
ature range the thermal diffusivity and thermal conduc- 
tivity are nearly proportional. At atmospheric pressure, 
thermal conductivity decreases with temperature, from 
the boiling point to the lowest temperatures at which 
thermal conductivity has been measured [22j. Water's 
behavior in this regard is unique among non-metallic liq- 
uids of low molecular weight, as all other such liquids 
show an increase in thermal conductivity upon cooling 
& 

Benchikh et al. have noted a strong correlation be- 
tween the thermal conductivity of both supercooled and 
stable liquid water at low temperatures and the ther- 
modynamic speed of sound c (in the limit of zero fre- 
quency) [2l|. We find excellent agreement between the 
experimental data and two classical formulations for the 
thermal conductivity of liquids, as shown in Fig. 3. The 
first is Bridgman's formula [25| as adapted for polyatomic 
molecules: 
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where k B is Boltzmann's constant and v is the molecular 
volume of the liquid, that is, the molecular mass divided 
by the mass density p. The second formulation is due to 
Eyring, with a correction from Eucken [26| : 



A = 2.8k B v 
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where 7 is the ratio of the isobaric heat capacity to the 
isochoric heat capacity, c p /c v . The speed of sound can 
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FIG. 1. Thermal-diffusivity data of Benchikh et al. [211 ] and 
Taschin et al. [20(, compared with the thermal diffusivity cal- 
culated from Bridgman's (4) and Eyring's (5) formulae. We 
use the TSEOS to evaluate the thermodynamic properties in 
these formulae [l^ . The inset shows simulation results of 
Kumar and Stanley for the TIP5P model [12]. 
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FIG. 2. Heat capacity data of Archer and Carter [l| and 
Angell et al. [24| . The solid curve shows the prediction of the 

TSEOS Fl3. 



be related to the adiabatic compressibility k s and the 
isothermal compressibility k t as follows: 
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allowing us to re-write Bridgman's formula in terms of 
thermodynamic properties as 



A = 2.8k n v~ 2/3 



1/2 



(4) 
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and Eyring's formula as 

/ , xi/2 

A = 2.8fc B ^ 2/3 — , (5) 
\pn T J 

We have used the TSEOS [ljj to evaluate these expres- 
sions. 

Bridgman's and Eyring's formulae differ little in the 
region where experimental data are available. According 
to two-state thermodynamics, both the isothermal and 
the adiabatic compressibilities show maxima associated 
with the existence of a virtual liquid-liquid critical point. 
These maxima are located close to the Widom line, de- 
fined as the locus of maximum fluctuations of the order 
parameter, a continuation of the liquid-liquid transition 
line into the one-phase region |14l - fl7| . We note that the 
magnitudes of these maxima are strongly dependent on 
the value of the critical pressure, the value of which (pos- 
sibly ranging from 10 to 40 MPa) is difficult to determine. 
In conclusion, the observed behavior of thermal conduc- 
tivity is anomalous inasmuch as it decreases upon cool- 
ing, it tracks the anomalous behavior of the adiabatic 
and isothermal compressibilities as shown in Fig. 3. 
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FIG. 3. Thermal conductivity calculated from thermal dif- 
fusivity [HHH and the TSEOS [13 ■ The formulations due 
to Bridgman (4) and Eyring (5) are both shown, as is the 
IAPWS correlation for thermal conductivity (which is only 
guaranteed down to 273 K) [z| • The inset shows the simu- 
lation data of Kumar and Stanley for the TIP5P model [12] , 
and the dashed curve on the inset is a quadratic fit to the sim- 
ulation data. Tw refers to the Widom temperature in both 
the main graph and the inset. 

Recently, Kumar and Stanley [l2| reported evidence 
of a thermal conductivity minimum in a simulation of 
the TIP5P model of water. The simulation results indi- 
cated that the thermal conductivity of this model at first 
decreases upon cooling as has been observed experimen- 
tally in real water; it then reaches a minimum at approxi- 
mately 255 K, and increases as the temperature is further 
decreased [12| . The latter increase might seem at first 



to conflict with experimental data for real water. How- 
ever, Kumar and Stanley's simulation locates the Widom 
temperature for the TIP5P model at roughly 245 K at 
atmospheric pressure, while thermodynamic equations of 
state place the Widom temperature for real water close 
to 228 K [l4T - tFi| . Rescaling the temperatures in the sim- 
ulation results so that the Widom temperature occurs at 
228 K places the predicted minimum at 237 K, several 
degrees below the lowest-temperature measurements of 
the thermal conductivity, so there is no real contradic- 
tion between simulation data and the predicted behavior 
of thermal conductivity in real water. Moreover, one can 
expect the minimum of thermal conductivity observed 
in this model to be smoothed by finite-size effects, as is 
typical for simulations. As can be seen in Fig. 3, both 
the Bridgman formula (4) and the Eyring formula (5) 
indicate that thermal conductivty should go through a 
minimum. We shall return to this topic in more detail in 
Sec. IV. 

Next, we address the possible effects of critical fluctu- 
ations on the thermal conductivity of supercooled water. 
It is well documented that the thermal conductivity of 
water diverges at its vapor-liquid critical point, and the 
associated anomaly affects the thermal conductivity no- 
ticeably throughout the critical region [27j ■ We investi- 
gate the possibility of such a divergence of the thermal 
conductivity near the liquid- liquid critical point of H 2 0, 
and any effects that this might have on the measurable 
behavior of the thermal conductivity in supercooled wa- 
ter. 



III. PREDICTIONS OF MODE-COUPLING 
THEORY 

In the vicinity of a critical point, couplings among the 
various hydrodynamic modes of a system become increas- 
ingly significant as fluctuations in the system become 
long ranged. This leads to anomalous behavior of the 
thermal conductivity in the critical region, including a di- 
vergence of the thermal conductivity at the vapor-liquid 
critical point. This divergence has been observed in strik- 
ing agreement with the mode-coupling theory in many 
molecular fluids near their respective vapor-liquid criti- 
cal points [13, HH . (This mode-coupling theory describes 
dynamics in the vicinity of a critical point and should 
not be confused with the mode-coupling theory of the 
glass transition). Due to a mode-coupling contribution 
to the thermal diffusivity, which arises in molecular fluids 
from a coupling between the heat mode and the viscous 
mode, thermal conductivity can be expected to diverge 
near any critical point at which the isobaric heat cap acity 
diverges more strongly than the correlation length £ [28| . 
Such a strong divergence of the isobaric heat capacity 
is a feature of the TSEOS [17], as well as other related 
scaling models EMI; this prompted us to investigate 
the possibility that a critical enhancement to the thermal 
conductivity might be experimentally observable. 
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In order to carry out our mode-coupling calculations 
we once again made use of the TSEOS [l7| for thermo- 
dynamic properties. This formulation is renormalized by 
critical fluctuations, and asymptotically close to the crit- 
ical point it is identical to the scaling models referred 
to above fl4T - [l6l |. For the background value of the ther- 
mal conductivity we used the formulation of the Interna- 
tional Association for the Properties of Water and Steam 
(IAPWS) [12]. IAPWS provides a formulation for the 
viscosity of water at atmospheric pressure that is valid 
for temperatures as low as 253.15 K [2!| |3(J. Extrap- 
olating below that value is a more subtle task. Data 
taken by Osipov et al. HH from 238 K to 273 K show a 
clear super- Arrhenius dependence on temperature ( "frag- 
ile" behavior in Angell's nomenclature [32|). Some re- 
searchers [33j have found evidence of Arrhenius temper- 
ature dependence ("strong" behavior) close to the glass 
transition, and thus of a fragile-to-strong crossover in wa- 
ter; such a crossover has been observed to occur at 228 K 
in confined water |34| . Other measurements, however, 
have found super- Arrhenius behavior at the glass transi- 
tion [35|]. Starr et al. have used Adam-Gibbs theory to 
estimate the viscosity in the experimentally inaccessible 
region 3(i . and this extrapolation includes a fragile- to- 
strong crossover as well. We have fit a super-Arrhenius 
equation to the experimental data of Osipov et al. (see 
Fig. 4), and in making the fit we have chosen a hypo- 
thetical temperature of structural arrest so that our ex- 
trapolation agrees with that of Starr et al. in the portion 
of the unstable region for which we make predictions. In 
our calculations we use the IAPWS formulation for vis- 
cosity at atmospheric pressure [2!| |3(| for temperatures 
above 254 K and our extrapolation for temperatures be- 
low 254 K. 

Mode-coupling theory gives a pair of coupled equations 
for the critical enhancements to viscosity and thermal dif- 
fusivity, and these equations should in principle be solved 
by iteration [37]. However, the viscosity anomaly associ- 
ated with critical fluctuations is very weak, so we work 
only to one-loop order in the iteration and simply use the 
background value of viscosity, which is strongly tempera- 
ture dependent (see Fig. 4). The anomaly of the thermal 
diffusivity is additive in nature, meaning that the thermal 
diffusivity can be split into a background and a critical 
contribution 12711: 



a b + Aa, 



(6) 



and the thermal conductivity can be treated similarly: 



A = pc p a = Ab + AA. 



(7) 



Even to one-loop order the integral for the thermal dif- 
fusivity enhancement cannot be solved exactly except in 
the asymptotic limit £ — ¥ 00. Our approximation strat- 
egy for treating the crossover from critical to mean-field 
behavior is based on the model put forward by Olchowy 
and Sengers [H, [3^] . It yields an expression of the form 

AA = pc p Aa = P c p ^f (n - fio) , (8) 
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FIG. 4. Viscosity data taken by Osipov et al. at atmospheric 
pressure (dots) [3l]]. The solid curve shows the IAPWS for- 
mulation for viscosity at atmospheric pressure [2S|, [3(J , while 
the dashed curve is a fit to the data of a super-Arrhenius or 
Vogel-Fulcher-Tamman law: 77 = 0.0885 exp [220/ (T - 197)], 
with T in K. The inset shows the product of viscosity and 
thermal conductivity. The increase in the viscosity completely 
dominates the decrease in the thermal conductivity, so it is 
clear that the decrease in thermal conductivity is not a result 
of the increase in viscosity. 



where R D is a universal amplitude very close to unity 
(experiments by Burstyn et al. [|o| hud Rd = 1-02 ± 
0.06). In the limit £ — > 00, we have (n — n ) — > 1. If we 
adopt R D = 1, this expression tends to the well-known 
limit of a Stokes-Einstein law for thermal diffusivity in 
which the correlation length of the critical fluctuations 
replaces the hydrodynamic radius of Brownian particles: 



Aa 



67T7/£ 



(9) 



Due to the effects of long-time tails on the hydrodynamic 
modes, mode-coupling effects do not completely vanish 
far from criticality [38|, |4l| . These long-time effects are al- 
ready present in the background, and so the phenomeno- 
logical term tt is introduced to subtract these effects 
from our expression so that it represents only the criti- 
cal enhancement. Further details of the approximation 
scheme can be found in the appendix to this article. 

The path along atmospheric pressure is not the crit- 
ical isobar, and along this path properties may exhibit 
finite anomalies but they do not diverge. We find that 
at atmospheric pressure, the critical enhancement to the 
thermal conductivity would reach its maximum in the 
vicinity of the Widom temperature and at a value close 
to AA = 8 x lO" 6 Wm-'K" 1 (see Fig. 5). This effect is 
certainly too small to be measurable, either in the exper- 
imentally accessible regime or in the region of the phase 
diagram below the homogeneous nucleation line (predic- 
tions for this region appear in Fig. 5 as a dashed curve). 
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Figure 5 shows experimental data as well as the predic- 
tion from Eyring's equation (2(|. Any effect from the 
critical enhancement is far too small to be visible on such 
a graph; it is shown in the inset, magnified by a factor 
of a million. We note further that error bars on thermal- 
conductivity measurements in water are typically of the 
order of Wm^'K" 1 , several orders of magnitude 

larger than any possible effect induced by critical fluc- 
tuations at atmospheric pressure. Even at the critical 
pressure (for which the TSEOS uses 13.1 MPa [13), and 
even if one could somehow carry out measurements in 
"no man's land" , any critical enhancement would be un- 
detectable, as it would be confined to to small a range of 
temperatures. 
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FIG. 5. The main figure shows thermal conductivity in water, 
both from experimental data [2(| Hl|] and from the Eyring's 
equation (5) evaluated with the TSEOS [13]. The inset shows 
the critical enhancement in the same units, but magnified by 
a factor of one million (10 6 ). The curves change from solid to 
dashed at the temperature of homogeneous nucleation. The 
critical effect is several orders of magnitude smaller than the 
background and will be completely undetectable. 



IV. DISCUSSION 

While thermal transport near the vapor-liquid criti- 
cal point is dominated by a Stokcs-Einstein law for ther- 
mal diffusivity, near the hypothesized liquid-liquid criti- 
cal point thermal transport will continue to be governed 
by the thermodynamic properties. We can identify two 
immediate reasons for this striking difference between 
the two critical regions. First of all, mode-coupling the- 
ory predicts that the critical enhancement to the thermal 
conductivity will be inversely proportional to the viscos- 
ity. The viscosity of supercooled water increases dramat- 
ically as the temperature decreases: Osipov et al. [3l| 
found an increase of a factor of 10 between 273 K and 
238 K in their experiment. Our extrapolation predicts 



that the viscosity at the Widom temperature will be yet 
another order of magnitude larger than at 238 K. This 
means that the viscosity of water at the Widom tem- 
perature for the hypothesized second critical point is be- 
tween two and three orders of magnitude larger than at 
the liquid- vapor critical point, which greatly suppresses 
any critical enhancement to the thermal conductivity. 

The second physical reason why any measurable 
anomaly will be confined to such a tiny region of the 
phase diagram in supercooled water is that while temper- 
ature is (to a good approximation) the thermal field for 
the vapor-liquid phase transition, for the hypothesized 
liquid-liquid phase transition in water it very nearly plays 
the role of the ordering field Near the critical 

point, the ordering field is related to the order parameter 
according to a power law <j> ~ foi^i with 6 « 4.8 for the 
Ising- model universality class [42|, [43|- Thus, asymptot- 
ically close to the critical temperature, small variations 
in the temperature correspond to large variations in the 
order parameter, and so for practical purposes a small de- 
viation from the Widom temperature moves the system 
far from criticality. 

The thermal-conductivity minimum merits some fur- 
ther discussion. We have noted, following Benchikh et 
al. [ll|, that thermal conductivity in supercooled wa- 
ter follows the thermodynamic (low-frequency) speed of 
sound in good agreement with Bridgman (4) and Eyring's 
(5) formulae for thermal conductivity in a liquid. Trinh 
et al. have raised the possibility of a minimum in the 
speed of sound [3], and calculations from the TSEOS 
predict such a minimum at 233 K [l7j | . So it is plausi- 
ble that thermal conductivity does indeed have a mini- 
mum, and that it continues to follow the speed of sound 
in the medium as it does throughout the measured part 
of the supercooled regime. The anomalous behavior of 
the speed of sound probably results from the existence of 
a region of the phase diagram in which it is thermody- 
namically rather inexpensive to convert water between 
the HDL to LDL structures, leading to a higher com- 
pressibility. Once this condition no longer obtains, more 
typical liquid behavior might be expected to resume. In 
explaining their simulation results, Kumar and Stanley 
have noted that many other systems show an increase in 
thermal conductivity upon ordering [l2j . It is interesting 
to note that strong links between thermodynamics and 
dynamics are a general feature of the supercooled region 
of glass-forming liquids; a good discussion of the link can 
be found in }45| . However, the correlation between com- 
pressibility and thermal conductivity in water is robust 
in much of the stable liquid region as well and it is not 
associated with glassy dynamics. 

Finally, we note the possibility of an even more fun- 
damental difference between water near its vapor-liquid 
critical point and water in the supercooled regime, one 
that points to a possible direction for further research. In 
the vicinity of the vapor-liquid critical point, the fluctu- 
ations in the density are not associated with the conver- 
sion of molecules between distinct molecular structures. 
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No or little structural relaxation takes place, and the 
decay of density fluctuations is completely governed by 
diffusion dynamics. In supercooled water, contrarily, the 
existence of distinct molecular structures seems likely to 
introduce into the dynamics a relaxation time for conver- 
sion between these structures 19]. Although liquid water 
at atmospheric pressure is non-dispersive down to tem- 
peratures of —15 °C, high-frequency sound (w > 1 GHz) 
below that temperature shows a positive dispersion rela- 
tion. This phenomenon has been widely attributed to a 
structural relaxation process taking place in supercooled 
water [46U48l |; in particular, it may be related to the 
conversion between the kinds of hydrogen bonding char- 
acteristic of HDL and LDL water. We emphasize that 
while the dispersion of sound near the vapor-liquid crit- 
ical point is solely an effect of the relaxation of critical 
fluctuations, the dispersion of sound in supercooled water 
is most likely a visco-elastic phenomenon [lj| . The char- 
acteristic time of this relaxation probably does not de- 
pend on the wavenumber of any fluctuation of the density, 
and such a relaxation will be present even in the mean- 
field approximation. Near the liquid- liquid critical point, 
both kinds of effects — diffusive decay of density fluctua- 
tions and structural relaxation associated with hydrogen 
bonding — may be present. Coupling between a relax- 
ational visco-elastic mode and a diffusive decay mode has 
been observed and analyzed in polymers both near to a 
critical point and far from the critical region [49]. This 
analysis may provide a point of departure for further in- 
vestigation of transport in supercooled water. 



fusive decay, the relaxation of this order parameter would 
be associated with a characteristic time that is indepen- 
dent of wavenumber. This would have far-reaching im- 
plications for various dynamic properties of supercooled 
water. 
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Appendix: Details of the Mode Coupling 
Calculations 



We start from the mode-coupling integral as given in 
Ref. M: 



Aa 



k B T 
(2tt)3 



fc- 2 sin 2 (<9) 



c p (0)v(k)(l + a(k)p/ V (k))- 



(A.l) 



The Olchowy-Sengers approximation [41| reasons that 
close to the critical point the term ap/rj is negligible. 
As noted we use the background value of the viscosity 
in our approximation, so after carrying out the angular 
integrals the integral to evaluate is as follows: 



V. CONCLUSION 

At high temperatures, density fluctuations associated 
with the vapor-liquid critical point noticeably affect the 
thermal conductivity over a fairly broad temperature 
range. Our calculations suggest, however, that this will 
not be the case for the hypothesized liquid-liquid critical 
point in supercooled water. The effect of density fluctu- 
ations on the liquid-liquid critical point is too small to 
be measurable. The Stokes-Einstein law that describes 
thermal transport in the vapor-liquid critical region will 
not be applicable to the thermal diffusivity in the liquid- 
liquid critical region. On the other hand, the thermal 
conductivity and thermal diffusivity of supercooled water 
are strongly correlated with the anomalies of the ther- 
modynamic properties associated with the existence of 
a liquid-liquid transition. The minimum of the thermal 
conductivity, found in simulations by Kumar and Stanley 
should also exist in real water, being associated with 
the maximum of compressibility and with the minimum 
of the speed of sound. 

Moreover, supercooled water may differ even more fun- 
damentally from water in the vapor-liquid critical region 
due to the possible existence of a non-conserved order 
parameter associated with the structure of the water 
molecules [5C| . Rather than obeying the dynamics of dif- 



Aa 



3ir 2 T] J c p (0) 



(A.2) 



In the TSEOS the isobaric heat capacity can be ex- 
pressed in terms of the relevant scaling variables and a 
few parameters of the system as follows: 



c p = -A(l + 6AP)T(0+l) 



1 



A 2 (l + 6AP) 2 T 2 Xl , 

where A and b are parameters from the TSEOS with the 
values given in the supplement to ref. [l7| . is a back- 
ground term representing the dimensionless heat capac- 
ity of pure HDL. The strongest divergence in the isobaric 
heat capacity is in the strong scaling susceptibility Xi, so 
for our third approximation we separate the term con- 
taining Xi from the rest of the heat capacity, and treat 
the remaining terms (the sum of which we shall call A) 
as having no wave-number dependence. We then use the 
Ornstein-Zernike approximation for the wave-number de- 
pendence of Xi '■ 



Xi(0) 

1 + fc 2 c 2 : 



(A.4) 



where £ is the correlation length characterizing the fluc- 
tuations of the order parameter. 
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Henceforth c p without any explicit wave-number de- 
pendence will refer to the hydrodynamic value, c p (k — > 
0). With that notation, the integral that we must evalu- 
ate for the thermal diffusivity takes the form 



Aa 



dk 



A 



1 



.4 



1 + £ 2 k 2 



(A.5) 



which yields the modified Stokes-Einstein law 



Aa 



k R T 2 



6irrj!; it 



c — A A 
— arctan(g D £) H g D £ 



(A.6) 



For convenience we define 



n= 2 - 

7T 



c — A A 
— arctan(g D £) H g D £ 



(A.7) 



This result is still not entirely satisfactory because in the 
limit of vanishing correlation length, that is, far from 
criticality, it does not vanish. In fact: 



.. k B T 2 
lim - — 

5->o 67rry4 ir 



A A 
— arctan(q D £) H g D £ 



k B Tq r 



The physical reason for this as summarized by Olchowy 
and Sengers in refs. [H, 5l| is that mode coupling is 
also responsible for the "long-time-tail effects on trans- 
port properties" t 38j . These effects are not critical effects 
and will be observed in the background, so if we want 
to find the effects due to critical fluctuations we should 
subtract off this remnant. For this reason we subtract 
the following term from tt: 



1 — exp 



1 + (1 - A/c p ) (qZ) 4 



(A.9) 



This phenomenological expression has the following lim- 
iting behavior: 



so that we have 



lim fl(x) = — , 
lim n(x) = 0, 



lim Aa = 0. 

5^0 



(A.10) 
(All) 

(A.12) 



The complete expression, using the above definitions, 



is 



Aa 



k B T 



(fi-Oo). 



(A.13) 



Because we are interested only in the critical enhance- 
ment to the thermal conductivity, for the correlation 
length we should use only the critical enhancement to the 
correlation length. We estimate the "background" corre- 
lation length £ b by observing the correlation length at a 
temperature T rcf = 315 K far from any critical point and 
assuming that the background correlation length in the 
system is proportional to the space between molecules. 
So we have 



(T rcf ) 1/3 

For the correlation length in our calculations we use 



(A.14) 



(A.15) 



where £ is the correlation length predicted by the TSEOS 
0. 

For the wave- number cutoff g D , we used a correlation 
identified by Perkins et al. between the wave-number 
cutoff and the amplitude £ of the correlation length 
anomaly [39j : 



On 1 = 3.683£ - 1.336 x lCT 10 m. 



(A.16) 
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